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I. 



INTRO DUCT ION- THE WHISTLER NO ZZL E 



The phenomenon of entrainment is a process through which a 
moving fluid increases its mass flow, by "grabbing up" the 
medium it is flowing in. The moving fluid is called the 
working fluid, and might be the primary water jet in an 
ejector pump. The medium might be any other fluid, say oil, 
in the case of the ejector pump. 

The entrainment process has found practical use for 
quite a long time. Analytical descriptions of the 
entrainment jrocess exist; however, none is general enough to 
describe all cases of practical interest. Even the 
relatively simple case of a steady jet in still air is 
rather difficult because of the many configurations which 
one would like to descri.be. Thus interest and research in 
this area continues, perhaps with more intensity than 
previously due to the recent interest in new V/STOL 
configurations which would employ this process. One such is 
the Navy's XFV - 12A, which is expected to derive a large 
amount of thrust through the augmentation, or entrainment, 
of its primary thrust using wing mounted ejectors. 

The success of concepts such as the augmenter wing 
depends cn many factors, an important one being the rate of 
entrainment of the primary jet, and another being its 
efficiency. This rate of entrainment is the rate of 
increase of jet mass flow with axial distance away from the 
jet nozzle. The rate of entrainment and nozzle efficiency 
are characteristics which vary with nozzle configuration. 

As mentioned earlier , there is a large variety of 
steady nozzle configurations. Lobe, slot, and hypermixing 



nozzles are some which have drawn recent interest for their 
entrainment capabilities. One should not neglect unsteady 
nozzle configurations, however, as such effects as swirling 
the primary flow, and oscillating the primary flow in a bi- 
stable manner (fluidic oscillator), also show promise in 
their ability to entrain. 

One of the newest nozzle configurations to surface is an 
unsteady nozzle. Perhaps the simplest in construction of 
all unsteady nozzles, it is illustrated in figure 1. The 
whistler nozzle, as it is called , produces a loud pure tone 
while in operation. Additionally, it creates a relatively 
large rate of entrainment. It is the whistler nozzle which 
is the subject of this thesis. 
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WHISTLES NOZZLE CROSS SECTION 



II . 



BE LAT ED PAST WORK 

The earliest work done on the whistler nozzle was by Hill 
and Greene (reference 1). In 1374, they cited the discovery 
of a new phenomenon, characteristic of a new device which 
they named the whistler nozzle. Their basic configuration, 
as illustrated in figure 1, consisted of a convergent 
section, followed by a section of constant area followed b 
a step change in cross sectional area to another section o 
constant area. The length of the last section could be 
varied tc produce loud pure tones of varying frequency, 
through the excitation of a standing acoustic wave in the 
final section, much as in an organ pipe. Hill and Greene 
observed that the whistler produced a greatly increased 
mixing rate, and decided that this was due to acoustic 
stimulation of the jet. Increased mixing rate c£ an 

axisymmetric jet due to acoustic stimulation is a process 
that had been observed earlier by Crow and 
Cha mpagne (ref erence 5). 

In addition to the basic axisy mmetric configuration, 
many other configurations were found to produce the 
whistling phenomenon. A few are illustrated in figure 2. 
Common to all of these was that the step in cross section 
had to extend completely around the jet. Also noted was me 
apparent presence of a separation reattachment cycle, 
occurring at the nozzle lip. 

Nc further study seems to have been performed concerning 
the whistler nozzle, save a report in 1975 by Hill and 
Jenkins (reference 2) , in which they reestablished the 
operating characteristics of the nozzle, and its ability to 
entrain more air. 
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Figure 2 - VARIOUS WHISTLFP CONFIGURATIONS 
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Hill and Greene’s work, has been summarized, since their 
original report, in comprehensive reviews of the ejector 
state of the art by Sc hum (reference 3 ) , and by Viets 
(reference 4). These reports call attention to the whistler 
nozzle's promise in rate of entrainment, and the lack of 
data regarding the whistler thrust efficiency. 

Although no further research directly involving the 
whistler is available, much work, has been done in the area 
of free jets, and entrainment, which might prove valuable in 
a study of the whistler nozzle. 



An obviously important area of related work is the 
sensitivity of free jets to acoustic stimulation. Prcbanly 
the best known and most definitive work in this area was 
performed by Crow and Champagne (reference 5) . By mounting 
a speaker inside the plenum chamber they managed to excite a 
free axis ym metric jet, at various frequencies. The effect 
of the sound was to produce a velocity fluctuation in tne 
nozzle exit plane. It was observed that the turbulent 
structure of the jet was sensitive to the acoustic 
excitation. Through the use of high-speed schlieren 
photography it was observed that under acoustic stimulation 
the turbulent structure of the jet developed a large-scale 
vortex structure, as contrasted with the small scale vortex 
sheet surrounding a normal jet. Attendant to the formation 
of the large-scale vortices, an increase in the mixing rate 
was also observed. 



Others, such as Bevilaqua and Lykoudis (reference 6) , 
have observed the turbulent structure in steady, free jets. 
The presence of a large - scale vortex structure has been 
observed in such jets at low Reynolds numbers. It appears 
that the largo — scale structure dissappears as Reynolds 
number increases. However, according to the findings of 
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Crow and Champagne, it may be reexcited through accusric 
stimulation. It may be conjectured from experimental 
evidence that the entrainment process is a combination of a 
large-scale "scooping up" of the ambient fluid, and a smaU- 
scaie nibbling away at the ambient fluid. These entrainment 
processes would correspond to the large and small scale 
vertex structures around a free jet. 

A further area of related research, perhaps not obvious, 
is the well known Coanda effect. The basic whistler 
configuration sets up a flow situation somewhat akin tc that 
of a bistable fluidic element, or fluidic oscillator. The 
difference between the two cases is that the fluidic 
oscillator csci.3J.ates between two positions, while the 
whistler must oscillate between many circumferential 
positions. That such an oscillation takes place seems 
indicated by the nozzle configuration and the separation 
reattachment cycle noted by Hill and Greene. 

Reference 10 provides a good summary of the mechanism of 
tne Coanda effect, as well as its characteristics regarding 
reattachment and sensitivity to sound. Though the following 
comments relate specifically to two- dimensional flows, it 
is thought that the general behavior and mechanism should 
carry over to a three- dimensional situation. 

Figure 3 illustrates the Coanda mechanism in 
two-dimensional flow. As shown, a jet is issuing fren an 
orifice and flowing between two adjacent walls. Because the 
jet will at some 'time be nearer to one wall than the ether, 
and because the jet is entraining air, a' pressure difference 
across the jet will develop that pulls the jet to the closer 
wall. This attachment forms a region of circulating flow. 
The point of jet attachment is determined by this region, in 
which a balance is struck between flow entrained by the jet 
and flow reinjected by the jet. 
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Figure 3 - COANDA EFFFCT 
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Anything that disrupts this balance causes the attachment 
point of the jet to shift. 

Observation has shown that the attachment point in 
Ccanda effect is sensitive to jet velocity and sound, within 
a range of Reynolds numbers. In the range of Reynolds 
number from 200 to 3000, under the effect of increasing jet 
velocity, or a sound input of increasing amplitude, the 
attachment point has been observed to move upstream. This 
upstream movement appears to approach a limiting position, 
beyond which the point will not move. 

Reference 10 relates the effects of jet velocity and 
sound on the attachment point to the point of transition 
from laminar to turbulent flow of the jet. Both increasing 
jet velocity and a sound input will move the transition 
point closer to the jet exit. This movement will increase 
the length of turbulent flow across the recirculation 
region. A turbulent jet entrains more air than a laminar 
jet, and thus the balance of flow in the recirculation 
region is disrupted, with a resulting upstream move of the 
attachment point. 
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III. 



SCOPS OF RESEARCH 



The intent of this research ms to provide a clearer 
insight into the workings of the whistler nozzle, and to 
develop an improved understanding of the ways in which one 
may effect control over the characteristics of a free jet. 

Specific areas of research were the nozzle thrust 
efficiency, the nozzle entrainment rate, and visualization 
of the lip interaction process described by Kill and Greene. 
Accordingly, experiments were conducted to obtain relevant 
data. Thrusts were measured on a thrust bed, mass 
entrainment measured in a device patterned after that 
developed by Kicou and Spalding (reference 7) , and high 
speed motion pictures taken of tufts about the nozzle lip. 
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Figure 4 - EHISTLER NOZZLE AND PLENUM 
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V . 



DESCRIPTION OF EXPERIMENTS 



A. THRUST EFFICIENCIES 



Thrust efficiencies were obtained by measuring nozzle 
thrust. A comparison was made between measured and 
calculated thrust. 

Test equipment consisted of a thrust bed and the nozzle. 
The thrust bed was configured to measure horizontal thrust, 
and slid on two tracks, mounted to a solid foundation. The 
nozzle plenum was clamped to the thrust bed which was in 
turn secured to the foundation through a lead ring, mounted 
with a load cell, the device that actually measured thrust. 
Air was supplied to the nozzle from a large tank, connected 
to the plenum through a regulator and a flexible tube which 
led vertically away from the plenum, so as not to bias the 
thrust measurements in any way. Calibration was performed 
with the use of a tray, connected to the bed through a cable 
and pulley arrangement, upon which one placed weights of 
known magnitude. Peripheral gear required included a 
voltage source, digital voltage meter, and Wheatstone 
bridge. Figure 5 illustrates the set up. Plenum pressures 
were measured using a mercury manometer. 
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Figure 6 - THRUST BED C A LIEF ATI OK CURVE 
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Calibration of the thrust bed indicated a certain amount 
of hysteresis in the system, as u\ay be seen in figure 6. 
Not serious at thrusts of about 6 pounds, it did introduce 
problems at low thrust levels, and precluded reliable thrust 
measurement at low Reynolds numbers. 

Thrust measurements were made for various collar 
positions, over a range of nozzle exit velocities. Recorded 
data included gage pressure in the nozzle plenum, thrust 
readings from the voltmeter, collar extension, and ambient 
pressure and temperature. 

Measured thrusts were plotted against pressure ratio 
across the nozzle. Figure 7 presents an example curve. 
Ideal thrust was calculated according to the procedures in 
appendix A. Thrust efficiencies were then obtained and 
plotted against nozzle pressure ratio, exit velocity, and 
Reynolds number (based on the nozzle diameter). 
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EXPERIMENTAL THRUST CURVE 



The principle cf operation is relatively simple. The 
entrained air of a free jet flows in at a right angle to the 
jet axis. The ent rainirient chamber simply blocks this inflow 
witn a circumferential membrane, and seals the ends with a 
plate at the nozzle and an orifice at the axial station 
where one wishes to measure the total mass flow. Ideally, 
all of the entrained air is supplied to the exterior of the 
membrane from a regulated source. As a free jet does not 
support any axial pressure gradient, the proper amount of 
air to supply through the membrane is indicated by the 
disappearance of any pressure gradient across tne orifice. 
One may regulate the membrane flow to obtain this point. 
These primary aspects of the device are illustrated in 
figure 8. 



Critical design aspects of the entrainment chamber, 
especially the orifice sizes, are treated in appendix fi. It 
is sufficient here to simply stare that one entrainment 
chamber was designed to perform measurement of the total jet 
mass flow at various axial locations. 



One primary difference between this chamber’s geometry 
and that of ethers (references 7 and 8) should be noted. As 
shown in figure 8, exit orifices were positioned in the 
porous cylinder at various distances from the nozzle. That 
portion cf the porous cylinder extending aoove the orifice 
plate was sealed by a fiberglass cylinder attached to the 
orifice. Thus, in this situation chat porrion of the jet 
beyond the orifice and inside the porous cylinder exit was 
not in a truly free-jet flow situation. This was as opposed 
to the reference 7 and 8 configurations where the exit 
orifice corresponded to the exit of the porous cylinder. 
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hoped that this chambers georaetry would not affect the 
chamber readings greatly, although some effect was to be 
expected . 

Peripheral equipment to the basic entrainment chamber 
included flew regulating, metering, and measuring devices. 
Figure 9 presents a photograph of the chamber and its 
peripheral gear. Separate air sources, both fed by a large 
tank, were supplied for the nozzle and membrane. Nozzle 
flow was metered through an orifice and controlled with a 
regulator. Eecause of the much greater mass flow required, 
membrane air was supplied through a venturi, followed again 
by a flow regulator. Pressures across the oritiee were 
measured with an inclined water manometer. Pressures across 
the venturi were measured with a mercury manometer. 
Pressure across the entrainment chamber orifices was 
measured with an aichohol micromanometer. Nozzle plenum 
pressure was measured with a water manometer. 
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Membrane air temperature was measured just after the 
venturi, using a mercury thermometer. 



Entrainment measurements were made at various axial 
locations for various collar positions. Recorded data 
included pressure difference across the venturi and orifice, 
pressure difference across the entrainment chamber orifice, 
nozzle plenum pressure, ambient pressure, air temperature in 
the membrane supply, and collar position. Standard 
procedure for eacn run was to vary the membrane airflow so 
that one obtained both positive and negative pressure 
gradients at the chamber orifice. This provided a more 
precise location of tne proper entrained flow rate. It is 
of interest to note that, due to hardware limitations and 
the large ancunt of air entrained by the jet at the furthest 
axial locations measured, nozzle exit velocities were 
limited to around 100 feet per second at the further axial 
locations. This also affected the number of useful 
locations for whistler measurements, as a strong whistler 
oscillation required a relatively large exit velocity, say 
300 feet per second. 
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C. FLOW VISUALIZATION 



Flow visualization of the nozzle jet lip interaction was 
performed in two ways: high-speed motion photography of a 
ring of tufts at the nozzle lip, for the purpose of 
discerning a pattern to the interactions; and a mixture of 
lamp black and kerosene inside the nozzle collar in order to 
observe the behavior of the separation reattachment 
location . 

Photographic equipment consisted of a Honeywell Fentax 
single reflex camera, a Eolex Supreme handheld 16mm movie 
camera, and a aedlake Hycam 16mm high-.speed movie camera. 
The tuft ring was simply constructed of cotton threads held 
by a masking tape ring, and a photo of the ring is presented 
in figure 11. Observation in the lampblack experiment was 
performed with the naked eye and no peripheral equipment was 
necessary. 

Photographic tests were to be done in three phases. 
High-speed still shots were first taken to determine the 
usefulness of the tufts for visualization, and to obtain the 
approximate required film speed for motion photography. 
Next, medium speed motion pictures, 64 frames per second, 
were taken to confirm the necessity of the high-speed 
camera. Finally, high-speed motion pictures were taken to 
visually control the tuft motion such that useful 
observations could be made. No particular collar settings 
or exit velocities were used in photographic work. The only 
requirement was to produce the oscillation phenomenon, which 
was never difficult. As it turned out, a lack of time 
prevented the completion of the high-speed photography 
plia se . 
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LAMPBLACK POOLING ILLUSTRATION 



VI. 



RESULTS and disc ussi on 



Nozzle thrust efficiency results are presented in figure 
13, plotted for various collar positions. The whistler 
oscillation was produced at each collar setting for the 
whole range of Reynolds numbers shown, although sudden 
frequency jumps did occur in this range. These frequency 
jumps did net appear to change the nozzle thrust or thrust 
efficiency in any discontinuous manner. Instead the trend 
was as indicated, progressively decreasing nozzle efficiency 
with increasing collar extension. The zero collar extension 
efficiencies were found to oe about 6 percent lower than 
those found in a previous report (reference 8) . This is 
believed due to the extremely simple, unflared configuration 
of the nozzle, which may be seen in figure 4. Reference 8 
nozzles were divergent nozzles and as such should be more 
efficient. Important to note is that the whistler 
oscillations may be obtained 'without much reduction in the 
basic nozzle efficiency. 
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from the chambers of reference 7 anrl 8. Additionally , in 
the case cf the whistler nozzle, it is believed the method 
of enclosing the jet up to the axial position of measurement 
may interfere with the sound - jet interaction (especially 
the cloth porous cylinder) and the formation of a large- 
scale turbulent structure. Figure 15 illustrates that 
readings taken were relatively constant with Reynolds 
number. As noted by Ricou and Spalding, the entrainment 
ratio can be a strong function of Reynolds number in certain 
ranges of Reynolds number. It is interesting to note that 
this behavior is supported by the present notion of jet 
turbulence structure and its development as exit velocity 
increases, which postulates that the major contributor to 
entrainment, the large scale structure, dissappears as 
Reynolds nu inter increases. 
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THRUST EFFICIENCY RESULTS 




Figure 14 - ENTRAINMENT RATIO RESULTS AND COMPARISON DATA 
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Figure 15 - ENTRAIN !1 ENT PATIO VS REYNOLDS NUMBER 
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LacK of time precluded testing at nore jet axial 
locations, as each location required a different orifice 
cylinder (see appendix A) . Of the original four orifice 
cylinders constructed using reference 3 data, only two were 
useable as the spreading angles encountered averaged about 
28 degrees, cr 1 0 degrees less than those used in reference 
8* In order to be effective, the orifices had to be 
positioned at axial distances further than they were 
designed for. One of those was modified to obtain readings 
at its design position, also . The larger two were never 
effective. Another problem encountered was the mass flow 
limitation cf the porous cylinder air supply. This 
restricted the useable nozzle flow at axial locations beyond 
X/D of about 6 to such low values that a whistler 
oscillation strong enough to affect the entrainment could 
not be produced. 

It was hoped that high speed notion pictures would be 
taken of the tuft experiment; however, this phase of the 
photographic experiments was not reached. The ef f ecti veness 
of the tufts in indicating the circumferential collar 
positions of attached and detached flow was proven in tne 
hign-spjeed still picture phase. Notion pictures with a fc>4 
frame per second camera showed that a much higher speed 
camera was necessary to slow the notion such that motion 
pictures would be useful in slowing the tuft motion to a 
reasonable speed. The necessary film speed was estimated as 
about 800 frames per second; however, the pictures were never 
taken. Seme representative still pictures are presented in 
figure 16. 

Lampblack tests appeared to be productive as indicators 
of the general behavior of the collar attachment point. 
Indications are that a sort of Coanda effect is involved, 
modified by the fact that while the jet is attached at any 
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point on the collar there are other points where detached 
flow allows a backflow into the collar jet yap, which tends 
to relieve the Coanda effect causing a circumferential 
shift in the point (or points) of attachment. At a fixed 
collar extension, with increasing jet velocity, the attachment 
point was observed to move upstream to a limiting position, 
just as the Coanda jet described in reference 10. However, 
at fixed exit velocity and increasing collar extension a 
different behavior was observed. The attachment point 
seemed to move with the collar until a point was reached at 
which the frequency of oscillation jumped. At that instant 
the point of attachment also appeared to jump suddenly back 
to its original position relative to the jet exit. This 
behavior indicated some interdependence cf the vjhistler 
attachment phenomenon and the oscillation phenomenon. 

In addition to the above experiments, two other 
observations were made which verified those of Hill and 
Greene. 

First, it was noted that extension of the whistler 
collar seemed to increase mass flow through the nozzle. 
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Figure 16 - TUFT PICTURES, RESULTS 
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This was confirmed with static pressure measurements in the 
nozzle 'throat 1 . These indicated that tne collar was acting 
as an overexpanded diffuser, resulting in an effective exit 
area slightly greater than the 'throat' cross section cut 
less than the collar cross section, and further resulting in 
a greater mass flow through the nozzle. 



Second, it was noted that a correspondence between 
frequency of oscillation, jet exit velocity, and collar 
extension seemed to exist. It was found that by gradually 
extending the collar as plenum pressure was .increased, one 
could 'follow' a frequency of oscillation. Accordingly, 
collar extensions were measured and jet velocity calculated 
at fixed frequencies of oscillation, and it was found that 
these fixed frequencies corresponded to constant ratios of 
collar extension to jet velocity, thus confirming the organ 
pipe character of the whistler oscillation noted i>y Kill and 
Gre ene. 
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VII. 



CONCLUSIONS AND H3CONMSN CATIONS 



The thrust efficiency and mass entrainment, 
characteristics of the whistler nozzle have been measured, 
in the case cf the entrainment, using a technique heretofore 
not used on an unsteady jet. Further, observations have 
been made which shculd contribute to a further understanding 
of the whistler flew mechanism. 

Several directions of further work are indicated. 
First, it is recommended that a thrust measurement device 
using a pendulum arrangement, similar to that used in 
reference 15, be constructed for a more accurate 
determination of the whistler efficiency. A device of this 
type might avoid the hysteresis problems encountered with 
low thrusts cn a conventional thrust bed. It is also 
recommended tnat a whistler nozzle of smaller throat 
diameter be constructed to permit mass entrainment 
measurements , at further jet axial locations, in the 
entrainment chamber, at the high exit velocities required to 
produce a strong oscillation. This will also permit further 
study of the possibility of whistler measurements in the 
chamber ever a greater range of Reynolds number. Next it is 
recommended that further flow visualization experiments be 
conducted. Specifically, high speed motion pictures of the 
tufts used here, and high speed scnlieren photography cf the 
jet, should prove useful. Finally, it is recommended that 
other unsteady nozzles be tested in the chamber, as the poor 
results obtained here with the whistler are only preliminary 
and may not be indicative of the usefulness of the chamber 
with non-acoustic effects such as swirl. 
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VIII. 



IOEoduction- the ping wtng 



The present work concerns a hollow cylinder, axially 
aligned, in supersonic flow and undergoing small amplitude 
oscillat ions in angle of attack. The intent is to calculate 
the resulting pressure distr ibut5_ons and aerodynamic moments 
on what might be mere properly termed a ring wing. 



Previous investigations, i.e. references 12 and 13, have 
applied the method of characteristics towards the separate 
calculation of the inner and outer flowfields, as well as 
pressure distributions and generalized aerodynamic forces, 
for such a cylinder in supersonic flow whose walls are 
undergoing small amplitude sinusoidal oscillations, or panel 
flutter. These previous investigations resulted in computer 
programs to perform the necessary calculations. 



The primary thrust of the present work was the 
inodif icat ion of the previously created programs and then the 
combination of their results to calculate the unsteady 



pressure distribution, lift and moment 
undergoing angle cf attack oscillations, 
modified programs, a parametric study of 
con ducted . 



on the cylinder 
Finally, using tne 
these items was 
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IX. 



PROBLEM FORMULATION 



This development is closely based on several previous 
papers , i.e. references 12 and 13. 

Consider a circular cylindrical shell, in supersonic, 
inviscid, adiabatic flow, whose axis is aligned witn the 
f reestream, and which is undergoing small amplitude 
oscillations in angle of attack, about the Z axis, (figure 
17) 

The governing equation for this flow situation is the 
linearized unsteady potential equation, in cylindrical 
coordinates. 



( 1-M 2 ) $ + <J> + 1 $ + 

xx rr r 

r 



$ 
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- 2M <i>. 



- 1 
c : 



= 0 
tt 



The relevant boundary conditions are the flow tangency 
condition and Sommerf eld * s radiation condition, i.e., that 
disturbances will propagate away from their source. 
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Figure 17 - RING WING ILLUSTRATION 
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Alpha, Oscilation Amplitude 



One further condition is necessar}^ along the X-axis, which 
amounts to setting the pressure disturbance at the axis 
equal to zero. 

The linearized form of the flow tangency condition is: 



M = lh + U 3h ,0 .< x $ l 
r=R 3t °° 9x 



which must be satisfied on the cylinder inner and outer 
surface. 

As stated, the boundary condition at the axis is a zero 
pressure disturbance: 




0 



The linearized form of the pressure coefficient is as 
fellows: 



C 



P 



-2 $ - 2 $> = P - P, 

U 2 t Uco q TO 

oo 



After nondimensionaliza tion, using cylinder length and 
length divided by freestream velocity, the above equations 
become: 
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(1-H 2 ) 4, xx + 



rr 



+ 1 



1 

T ; 



*06 - 2m2 Kt - 



tt 



C p “ -2( + ) 



subject to: 



<|> | = + 3_h ,o ^ x ^ 1. 

r=R 3x 



One now assumes a cylinder oscillation of the form 



h(x,6,t) = Z(x) cos(6) e 



ikt 



where the non-dimensional frequency, k, is: 



k = u 1 

Uoo 



The deflection amplitude, Z (x) , for the case 
consideration, is of the form: 

Z (x) = a ( x - x ) 



un ler 
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Following the assumed form of the deflection, 
perturbation potential must be of the form: 

■j ■£* 

<j>(x,r,0,t) = (t(x,r) cos(G) e 



and, after substitution, the basic equations become: 



(1-M 2 ) + <J> + 1 4> + (k 2 M 2 -1 ) 4> - i 2kM <j)^ = 0 

xx rr r r 



C p = " 2{ ^x + ik( ^ ] 



subject to: 




Z + ikZ 
x 



,o $ x 4 1. 



and : 



C 



r=0 



0 



the 
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X. METHOD OF CHARACTERISTICS 



For supersonic flow, the linearized perturbation 
potential equation is hyperbolic, and has characteristics 
which satisfy the following differential equation: 



( M 2 - 1 ) d 2 r - d 2 x = 0 



Since the cnaracteristics in supersonic flow correspond 
to the i-iach lines, then for as, the differential arc length 
along a characteristic: 

dr = + _1_ , dx = y/M 2 -1 

ds M ds M 

An arbitrary function, F(x,r) , has the following 
derivatives along the characteristics: 

dF = F dx + F^ dr = /m 2 -1 F v + 1 F„ 

ds x r di FT ” m r 

Defining dsl and ds2 to be the differential arc lengths 
along the left, and right-running Mach lines respectively, 
then : 




F 2 = dF 
ds 



2 
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The cross derivatives become: 



-F 



12 



= -F 



21 



= 

M‘ 



( ( 1-M 2 ) F 



xx 



+ F ) 
rr 



Solving now for Fx and Fr: 



F = 
x 



M 



2 /M- 1 



(F. 



F 2> 



F = M (F, - F,) 
r 2 1 2 



Now the basic equations may be written in terras of d> 
and its derivatives along the characteristics as follows: 



<f> 



12 



Or 




= 4> = 1 (4>, “ 40 + (k 2 

21 2rM 1 2 

= M (4»_ - Cf, ) = 2 + ikz 

2 1 2 x 

= —2(<fi + ikcf> ) 

x 



1 ) <J> - ikM (<f> + <f> ) 

r 2 M 2 M 2 -l 1 2 
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finite difference 



These equations have been put into 
forra and programmed in Fortran IV for the separate cases of 
exterior flow and interior flow of a cylinder whose vails 
were undergoing small amplitude sinusoidal oscillations. In 
the interior flow program the X axis boundary condition is 
satisfied on a small inner cylinder in place of the axis, to 
circumvent problems encountered should r go to zero. For 
detailed accounts of the equations and programs refer to 
references 12 and 13. 

The primary work here was the transformation of the 
aforementioned programs to the case under consideration , 
i.e., a cylinder undergoing small amplitude oscillations in 
angle of attack, and then a parametric variation study of 
the pressure distributions and generalized aerodynamic 
force. The programs were modified such that all cases run 
were with an angle of attack amplitude of 0.1 radians. 
Since the basic equation is in a linearized form, the 
conversion of any results to other angle of attack 
amplitudes is a simple matter of multiplication. For those 
interested, the programs resulting from changes required to 
the reference 12 ana 13 programs are listed in appendix D. 
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XI. 
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As stated above, the results studied were the pressure 
coefficient amplitude distribution, and the generalized 
aerodynamic force, Qmr. In its most general form, Omr is as 
follows: 



Qmr = 1 SC 'i' (x) dx 
~T Pra 



Z (x) = Z d 1 . (x) 

j 3 3 



where Cpm is the pressure coefficient amplitude resulting 
from the ra'th axial deflection mode, PSIr is the r*th axial 
deflection node, and A is the amplitude of the j * th axial 
deflection mode. For th* special case under consideration, 
Q takes the f ollovri r.g form: 

Q = 1 /C ¥(x) dx 
2 p 



Z (x) = H' (x) = a (x - x Q ) 
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Thus , in the case under study Q has the special 
interpretaticn of being the moment coefficient amplitude 
about the point Xc multiplied by the angle of attack 
amplitude, alpha, which is equal to 0.1 radians in all cases 
presented . 
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XII. 



RESULTS and disc uss ion 



Several cases cf comparison were considered towards the 
evaluation cf program results. Figure 18 presents two cases 
of comparison for pressure distributions. The outer radius- 
to- length ratio was taken to such a value (10) as to 
approximate a flat plate locally. The inner 
radius- to-length ratio was taken to such a value (9.75) as 
to cause a reflection on the inner surface of the outer 
cylinder. 

For sufficiently large radius - to - length ratios the 
exterior pressure distribution should approximate the 
theoretical value for a flat plate in supersonic flow, i.e. 
2a/ 8- Figure 18A presents this comparison, and as can be 
seen agreement is good. 

For sufficiently large radius - to —length ratios the 
interior pressure distribution should approximate the 
theoretical value for a flat plate. Additionally, if tno 
inner cylinder is of such a radius-to -length ratio as to 
cause a il a c h wave to reflect upon the outer cylinder inner 
surface, a flow situation exists which should approximate a 
flat plate in a free jet, where the mach waves are reflected 
from the jet boundary. This is due to the boundary 
condition imposed on the inner cylinder, C'p = 0. Reference 
16 has treated this case in analytic form for small 
frequencies. Figure 18, A and B, presents a comparison of 
program results and the results of reference 16. As can be 
seen agreement is good. 
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Figure 18 - AXIAL PRESSURE DISTRIBUTION, 2D LIMITING CASE 
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Figure 19 - AXIAL PRESSURE DISTRIBUTION, COMPARISON WITH 

ZIEREP 
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Reference 14 has treated the case of a nonoscillating 
ring wing in supersonic flow. Figure 19 presents a 
comparison of procrram results and those of reference 14, for 
the case shown. A s can bo seen agreement is fair. 

In reference 12, which developed the basic inner flow 
field program used in this report, the occurrence of a 
numerical instability in the pressure distributions, under 
certain conditions, was noted. It was found that this 
disturbance, which occurred only after a reflection, was due 
to an undefined constant in the subroutine which computed 
the boundary condition at the axis. After defining the 
constant correctly, the pressure distributions were se«n to 
settle down quite nicely. Figure 20 presents a before and 
after comparison of pressure distributions for a 
representative case. 
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Figure 20 



NUMERICAL INSTABILITY ILLUSTRATION 



Figures 21 through 23 present various parametric studies 
of the moment coefficient, Q, for cases Loth with and 
without a reflection on the inner cylinder surface. In 
these figures the coefficient plotted is actually the 
difference between that calculated for the inner cylinder 
surface and that calculated for the outer cylinder surface, 
divided by the angle of attack amplitude. Thus the 
coefficients plotted are total moment and pressure 
coefficients, as opposed to those plotted earlier, which are 
the coefficients for the inner and outer cases separately. 
As can be seen from figure 22, for a center of rotation at 
the leading edge the cylinder has definite geometric 
regions of stability with respect to angle of attack. 
Figure 21 illustrates the effect of increased frequency on 
the moment coefficient. Figure 23 illustrates a typical 
pressure distribution including one reflection, from which 
one can infer the reason for the behavior of Q in figure 22, 
and also the behavior of Q with a varying center of 
rotation . 
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Figure 22 - PITCHING MO MSN T YS OUTER RADIUS-TO-LENGTH 

RATIO 
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R/L = . 3 
r /L= ,0001 
U M =1.5 
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Figure 23 - REPRESENTATIVE AXIAL PRESSURE DISTRIBUTION 

INCLUDING ONE REFLECTION 
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XIII. 



COLLUSIONS AND PSCONNENDATIONS 



The computer programs of references 12 and 13 have been 
successfully adapted to the case of a ring wina in 
supersonic flow undergoing oscillations in annle of attac v . 
Available cases for comparison have shov/n good agreement, 
indicating the validity of the method. The programs irate 
possible the relatively easv calculation of the unsteal v 
lifts and moments of a ring wing in supersonic flow, for 
such applications as tubular artillery shells (STD 0 , 
spinning tubular projectile), although spin may induce 
aerodynamic effects not accounted for in the method. 
Further evaluation and parametric studies are recommended 
for a more complete verification of the computer programs. 
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APPENDIX A 



TSENTROPTC THRUST CALCULATIONS 



Isentropic thrust calculations made use of standard 
isentrooic relations. Calculations were based on Dlenum 
pressure and nozzle exit area. The formula for isentropic 
thrust is arrived at as follows. Nozzle thrust may 
written in the following form: 

T = w V =(pAM 2 JTRT) 
s ex ex 

where W is nozzle mass flow, A is exit area, and one has 
applied the speed of sound relation. Exit Mach number may 
be written in terms of pressure ratio as follows: 



M 



ex 



2 ( (P n /P ) * V* -1) 

y-i ex 



Applying this and the state equation to the formula for 
Ts one obtains: 



-1 



T 

s 



P A 
ex ex 



2 

im 



Now, once one 
and thrust, Texp, 
efficiency: 



((P /P ) -1) 

u ex 

has measured plenum and ambient pressure, 
one can calculate the isentropic thrust 



U 



T 



/ T 

exp 



s 



67 



APPENDIX B 



ENTRAINMENT RATIO CALCULATIONS 



Mass flows and entrainment ratios were calculated 
according to the equations and procedures set out in 
reference 11, the ASME Power Test Codes. The equation f or 
mass flow, from reference 11 is: 

w = 359 c f d 2 Fa Y /hp , lb/hr 
where : 



c Discharge Coefficient 

f =l//l-6 !( 3=d/D 

d and D d El "] j D 

^ 1 l 

^a Thermal Expansion Factor 

y (Ya for a venturi) Adiabatic Expansion Factor 
h The pressure drop across the device, in inches of 

water 

p Density of the fluid being metered, pounds p a r 
cubic foot 



If wl is the entrained flow 
uses a venturi, w2 an orifice) : 

Ei = Si -£i ( ^i )2 E|l ~1 

w 2 C 2 ^2 ^2 Fa 2 Y 2 



and w2 the primary. 



/h 

"E 



1 

2 



then (wl 
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Non-geometric factors for the primary flow, such as 
discharge coefficient, could be reasonably approximate! for 
all tests as: 

C2 = . 65 Fa 2= 1 . 0 Y2=.99 



In the case of the entrained flovr these factors had to 
be changed for each test axial location, due to the widely 
varying flow requirements. The Geometric factors above had 
the following fixed values: 

F1 = 1 . 0 32 3 d1 = . 75 F?= 1 . 07 92 d2 = 1.11 



Now, entrainment ratio is commonly written as: 

M = W-+Wp s =w^ + 1 

M 0 w 2 -r 2 

Using the previous formula for w1/v2, the entrainment 

ratio may now be written as: 

M = A /h, +1 

M 0 TT^ 

where A is the combination of the ratios of factors 

previously listed. This equation was used to calculate the 
nozzle entrainment ratios. 
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ENTRMNK2KT CHAMBER DESIGN 



The entrainment chamber design v.-as closely based on 
previous wort by Ricou and Spalding (reference 7) r and 
Peschke (reference 3). Peschke's report especially was a 
good source of specific design data, while Ricou and 
Spalding was primarily a source of the design principles. 

The entrainment chamber consisted essentially a 
central porous cylinder, closed at one end with a solid wall 
through which the nozzle protruded, the other end being open 
to the outside. Orifice plates of specific diameters were 
constructed and mated to cylinders of specific length, 
length and diameter specified according to the data 
available in Peschke. These orifice cylinders were of such 
outside diameter that they fit closely inside the oorous 
cylinder, into which they were inserted for the measurement 
of mass flow at the various jet axial positions 
corresponding to the orifice diameters. I’he purpose of the 
orifice was to seal the open end such that entrained air 
could only enter through the porous cylinder. The purpose 
of the orifice cylinder was to seal off that portion of the 
porous cylinder beyond the jet axial position where flow was 
being measured. The porous cylinder was surrounded 
everywhere, except at the open end, by a box, also encasina 
the nozzle plenum, which served as the nenbrane air supply. 

The entrainment, chamber was about 4 by 4 feet square and 
4 feet high, supported on 3 foot legs. The porous cylinder 
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was 30 inches in diameter and 2 5 inches long. Orifice 
diameters used, corresponding to the listed iet axial 
positions, were as follows: 



Diameter 


X/D 


1 . 75 


2.5 


3.0 


5.0 


6 . 5 


13.0 



The entrainment chamber box, legs, orifice olates arid 
porous cylinder £ rame work were all constructed of wood. 
The orifice cylinders were made of fiberglass sheet and. the 
porous covering to the porous cylinder frame work was two 
layers of dacron crepe cloth, as suggested in Peschke's 
report. 

One critical chamber design parameter was orific Q 
diameter. T t is clear that too small an orifice would 
interfere with the jet being measured, while too large an 
orifice would not properly restrict the entrained flow to 
entering only through the membrane, thus preventing the 
creation of a readable pressure gradient across the orifice. 
The approach to this problem in both references 7 and 8 was 
to find an optimal orifice size a-*- a specific jet axial 
location by testing several orifice sizes. To minimize 
orifice construction, this investigation took a different 
approach. That was to use the reference 8 orifice data as a 
starting point, and vary the orifice axial location, such 
that a readable pressure gradient was created while not 
interfering with the jet. In both approaches the attempt is 
to optimize the orifice diameter, jet axial location 
combination . 
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Oris primary difference between this chamber's geometry 
and that, of others (references 7 and 8) should be noted. As 
shown in figure 8, exit orifices were positioned in the 
porous cylinder at various distances from the nox 7 le. That 
portion of the porous cylinder extending above the orifice 
plate was sealed by a fiberglass cylinder attached to the 
orifice. Thus in this situation, that portion of the dot 
beyond the orifice, and inside the porous cylinder exit, was 
not in a truly free jet flow situation. This was as opnosed 
to the reference 7 and R configurations where the exit 
or.irice corresponded to the exit of the porous cylinder. 
Thus in those experiments the jet past the orifice was in a 
truly free jet flow situation. The configuration used in 
this report. was adopted in order to facilitate the 
measurement of mass flows at many axial locations. it was 
hoped that this chamber's geometry would no f affec*: The 
chamber readings areatly, although some effect was to be 
expected . 

An important, problem, again with regard to the chamber 
orifices, concerned the measurement of the pressure 
differencial across the orifice. Roth reference 7 and 8 
indicated that a measurement accuracy of .001 inches ot 
water was required. Reference 7 used a micro manometer 
while reference 8 used a pressure transducer for this 
purpose. After the entrainment chamber was constructed 
experience showed that such accuracies were required. This 
investigation used an alchohol micro manometer whieh 
provided an accuracy of .0001 inches of nanometric fluid. 
This accuracy was obtained throuqh the combination of an 
inclined tube, hairline indicator, magnifying lens and a 
micrometer adjustment to indicator level. 

Flow measurement devices for the primary and entrained 
flow were an orifice and venturi respectively. These were 
constructed and used according to ASMS Power Test Co^es 
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specifications (reference 11). Based on these 
specifications, an incline} water manometer was f ounl to he 
necessarv for measuring the orifice pressure drop, while a 
vertical mercury manometer was sufficient for the venturi. 
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APPENDIX D 



PROGRAM LISTINGS AND DOCUMENTATION 



This appendix contains basic program and input 
descriptions for the two programs used, f rom references 12 
and 13. The programs also are listed in their modified form 
for the computation of a ring wing at angle of attach. The 
program statements relevant to the modification are 
indicated in the program listings. Also indicated in the 
interior flowfield program is the statement vfhose omission 
caused the numerical instability noted in reference 12. 

The inner flow field program consists of 8 subroutines, 
and a main program which manipulates each of these. The 
action of each particular subroutine is indicated in the 
program listing. Inputs are made by way of a Fortran 
Namelist. The input parameters are, in proper order: 

FSTEflN Free stream mach number. 

RFDFRQ The reduced frequency. 

RO The outer cylinder radius to length ratio. 

I 

RI The inner cylinder radius to length ratio. 

NGRDFN Grid fineness, or the number of grid points 
taken on the first right, running mach line. 

N Circumferential mode number. 

M Axial mode number. 

Ml The m in Qmn, the generalized aerodynamic force. 
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IPSTNT ?.n output parameter. 0 indicates only basic 
output. 1 indicates output including such things as Phi and 
its derivatives at each flow field point. 



Of the above 


inputs, in the 


modi f ied 


program , 


t he 


i nputs 


N, and Ml 


must be set 


to 1 , 


as the 


ca so 


under 


consideration is 


not panel flutter but a 


cylinder 


at 


angle 



of attach. The other inputs are arbitrary. 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

CYLINDRICAL SHELL/PANEL FLUTTER PROGRAM 

ccccccccccccccccccccccccccccccccccccccccccccc cccccccccccccc 

THIS IS THE MAIN PROGRAM 



COMPLEX PHI (400} 
DIMENSION X( 400 ! 



» D p HI S 1 ( 400 ) 5 D PH IS 2(400) 
t R (400 ) t X P T (400) ,CPR(400) ,CPI(400) 
DIMENSION DUMMY ( 400) t DUNCE (400 I 

C CM. V *0N/ BLCK1 / FSTRHN, DELTAS » OX, DH,DSTSTR ,R0 ,RI tREDFRQ 
A 1 , A 2 , A3 » A4 t R N * R M » R M 1 * F AC T 1 1 FACT2 
NGRDFN »NPTS» LCCUNT t IHAVEP t M ,M1 1 I PRINT 
PHI t DPH I SI » DP H I S 2 
X, R,XPT ,CPR,CPI 



COMMON/ SLCK2 / 
COMMON/ BLCK3/ 
COMMON/ BLCK4/ 
COMMON /BUCKS/ 



DATA PI/3.141593/ 

ENDPT ( Y 2 » Y 1 1 XO , X 2 i DELX J = 
1 Y2 



( Y2 - Y 1 )*( XO -X 2 ) /DELX 



EPS = 1 cE-04 
I ER = 0 

CALL INPUT! I ER) 

IF ( IER .EQ. i ) GO TO 12 
1000 CONTINUE 

IHAVEP = 1 
LCCUNT = 1 

CALL I N IT AL 

1 CONTINUE 
CALL COMPXR 

IS THE GRID POINT ON THE INITIAL RIGHT RUNNING MACH 
LINE? 

IF (( LCCUNT .EQ. 1 ) .AND. (IHAVEP .LE. NGRDFN )) 

1 GO TO 3 

IS THE GRID POINT ON THE OUTER CYLINDER SURFACE? 

IF ( IHAVEP .EQ. 0 i GO TO 5 

IS THE GRID POINT ON THE INNER CYLINDER SURFACE? 

IF ( IHAVEP .EQ. NGRDFN ) GO TO 6 

THEN THE GRID POINT IS IN THE GENERAL FLOW FIFLD. 

2 CALL GENFPT' IER) 

IF ( IER .EQ. 1 ) GO TO 12 
IHAVEP = IHAVEP * 1 
GO TO 1 

3 CALL MACHLN 

IF ( IHAVEP .EQ. NGRDFN ) GO TO 4 
IHAVEP = IHAVEP <- 1 
GO TO 1 

THIS IS THE LAST POINT ON THE INITIAL MACH LINE. 

4 LCCUNT = LCCUNT + l 
IHAVEP = 0 

GO TO 1 

5 CALL RAD2 

IS THIS GRID POINT PAST THE FLEXIBLE CYLINDER LENGTH? 
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IF ( ( X < 1 ) +EPS ) .GE. I. ) GO TO 8 
IHAVEP = IHAVEP f 1 
GO TO 1 
C 

6 CALL RADI 
IHAVEP = 0 
LCCUNT = LCOUNT + L 
GO TO I 
C 

8 CONTINUE 

IF ( ABS(XPT(LCOUNTi - 1. ) .LE. EPS ) GO TO 10 
L = LCOUNT 
LM1 = LCOUNT - I 
L M2 = LCOUNT - 2 
XFT(L) = l. 

C P R ( L } = ENDPT ( C PR (LM 1 i ,CPR(LM2) , 1. 00 , X P T( LK1 ) ,DSTSTP) 
CPKLJ = ENDPT (CPI (LM1) ,CPI (LM2 ) , 1 . 00, X PT ( L Ml ) , DST STR i 



WRITE(6» 105) 

00 9 I = I , L 

ARG = RM1*PI*XPT< I ) 

THE FOLLOWING TWO CARDS WERE MODIFIED FROM THE 
ORIGINAL REFERENCE 12 PROGRAM. 



DUMMY ( I )=CPR ( I )*. 1 *XPT( I ) 

9 DUNCE ( I }=CPI ( I)*.1*XPT(1) 

SUM I = 0 . 5=M CPR ( L ) «• CPR(LMl) )*< 1.00 - XPTU.MU) 
SUM2 = 0«5*<CPI(L) * CPI(LMl) )*( 1.00 - XPT(LMl)) 
C 

WR I TE(6?L14) FSTRMN, REDFRQ , RO , R I , N 
C 

CALL QSF ( DSTSTR , DUMMY » DUMMY » LM 1 ) 

CALL QSF (DSTSTR, DUNCE , DUNCE »LM1 l 
CRFAL = 0.5*<PUMMY(LMI) * SUM!) 

Q I MAG = 0 • 5* ( OUNC E < LM 1 ) *• SUM2) 

C 

WRITE ( 6, I 15) M, Mi , QRE AL , M » M 1 » QI MAG 
GO TO 12 
C 

10 CONTINUE 
XPT(L) = 1. 

WRITEC6* 105) 

DO 11 1 = 1, LCOUNT 
ARG = RM1*PI*XPT < I ) 

DUMMY ( I ) = CPP.(I ) *. 1*XPT( I) 

11 DUNCE ( I) = CPI( I )*.1*XPT( n 

c 

WRITEt 6 , 1 14) FSTRMN, REDFRQ, RO, RI ,N 



CALL QSF( DSTSTR, DUMMY, DUMMY, L) 

CALL QSF (DSTSTR, DUNCE, DUNCE ,L) 
QREAL = 0 * 5*DUMMY ( L i 
QI MAG= 0.5*DUNCE(L) 

WRITE (6, 115) M, M l , QREAL , M , M 1 ,Q I MAG 



12 



998 



999 



150 
105 
110 
1 14 



CONTINUE 
WR I TE ( 6 , 99 8 ) 

FORMAT ( • 1* ,’CP REAL' ) 

CALL PLOTP(XPT,CPR, LCOUNT, 0) 
WR ITE( 6 ,999) 

FCRKAT( • l ' , • CP IMAGINARY’) 
CALL PLOTPtXPT, CPI , LCOUNT, 0) 
STOP 

FCPMATt FI 0.5) 



,T19, ' I ' ,T27 , ’ X 1 



FORMAT ( 8 1’ . , . - 

FORMAT ( 8 0' ,15X,I5,3(3X,F9.5!) 
FORMAT ( // ' 

1 'REDUCED 

2 ’LENGTH) 



, T 3 8 , ' CPR’ , T50 , ' C PI ’ ,//) 



• , T 2 0 , 1 MACH NUMBER = 
FREQUENCY = ' , F10.7, 
OUTER = ' , F 1 0. 7 ,/• ' 



/ * 



,f 10 



7,//’ ' , T 20 



• , T 2 0 , ’ (RADIUS/' 



T20 , ' tRACI US/LENGTH) 



I 
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3 ,' INNER = *,F10.7,/' • , T20 , • C I RCUMF E RENT I AL MODES 

A /• • ,T34, • NUMBER = ',15,//) 

I I 5 FORMAT! 'O' ,T20, • Q{ • , I 2, • , « , 12, ‘ )REAL = ' ,F l 0. 5, / ' 0 ' , 

I T20, «Q( ' , 12,' , ' , 12, • ) IMAG = SFI0.5,/) 

136 FORMAT! » 1 ' ) 

END 

THIS IS SUBROUTINE INPUT 
SUBROUTINE INPUT! I ER) 

SUBROUTINE INPUT READS INPUT VARIABLES, AND DEFINES 
SEVERAL CONSTANTS USED THROUGHOUT THE PROGRAM. 

C CMMCN/ BLC K 1 / F S TR M.N , DE L TA S , D X , DH , DSTSTR , R 0 , R I , RE DF P . Q 
COMMON/ BLCK2/ A1 , A 2 , A3 , A4 , R N ? R ’ ! , RM I , F A C T 1 , F AC 7 2 
COMMON/ BLC K3 / NGRDFN , N PTSi LCCUNT , IHAVEP , N, M,M1 , I PRINT 
NAMEL I ST/ NAM I / F S TRMN , REDFRQ ,RO, RI , NGRDF N , N ,M , ML , 

1 I PRINT 

WRITF(6 ,80 ) 

READ! 5, KAMI ) 

WRITE (6.NAK1 ) 

FMANGL = ARSIN! L . /FSTRMN) 

DH = < R 0 - RI )/FLOAT( NGRDFN) 

IF ( DH . LT. 0. > GO TO 1 

DELTAS = DH*F S TR MN 
DX = DH /TAN! FMANGL) 

DSTSTR = 2.*DX 
NETS = 1. /DSTSTR + 1 

IF ( NPTS .GT. AGO ) GO TO 2 

BETA = SORT! FSTRMN^FSTRMN - 1.5 
RN = FLCAT(N) 

RM = FLOAT (M) 

RF1 = FLOAT (Mi) 

FACT 1 = REDFRQ^R EDFRO 

FACT2 = RN*RN/(FSTRMN*FSTRMN) 

COMPUTE CONSTANTS TO BE USED IN THE COMPUTATIONAL 
MOLECULES . 

A I = 0. 25*DELTAS/FSTRMN 
A2 = 0 • 5*F STRMN/ BETA 
A3 = REDFRC=~DEITAS*A2 
AA = 0.5-F STRMN 
RETURN 



I I FR = 1 

WRITE! 6 ,100) 
GO TO 3 



2 I ER = 1 
WRITE(6» 1 10) 

3 RETURN 

80 FORMAT! * 1* ) 

100 FORMAT!//' 35X ,' ABNORMAL 
1 ' INCORRECT' ,// ) 



TERMINATION: INPUT RADII ', 



110 FORMAT!//' 
1 *E THE 
END 



'• , 35X , • ABNORMAL TERMINATION: 
SIZE OF VECTORS XPT, CPR, AND 

THIS IS SUBROUTINE INITAL 



MUST 
CPI * 



INCREASE 

,//) 



SUBROUTINE INITAL 

SLBROUT INE INITAL INITIALIZES ALL FLOW FIELD 
QUANTITIES. 

COMPLEX DPHIDR, DPHIDX 
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CCMPL 5X PHI UOC)) , DPHI SI (400) ,DPHIS2(400) 

DIMENSION X( 4001 » R ( 400 ) ,XPT(400) ,CPR(400 ).CPI(400) 

CCMM0N/8LCK1/ FS TRMN, DE LT AS t DX , DH , DSTST R , R0 , R I , RE DFRQ 

CCMM0N/BLCK2/ AI , A2,A3,A4,RN,RM,RMl ,FACT1 ,FACT2 

CO MON/ BLCK3 / NGRDFN » NPTS » LC GUNT , I HAV EP , N , M , M 1 , I PR I NT 

COMMON/ 6LCK4/ PHI ,DPHI SI , DPHIS2 

COMMON / S LC K5 / X , R , XP T , CPR , C P I 

DATA PI/3.141593/ 



X ( 1 } = 0. 

R ( 1 > = RO 
X FT ( 1 1 = 0 . 



ALONG THE INITIAL RIGHT RUNNING CHARACTERISTIC PHI AND 
DPHIS2 ARE ZERO. 

PHI11) = (0.,0.) 

DP H I S 2 ( 1 ) = (0..0.) 



AA=. 2/F STRMN 
DPHI SI ( 1 J = C M P L X ( AA»0. ) 

DPHIOR = A4*DPHISl U> 

DPHI DX = A 2* DPHI SI i 1 ) 

CPR 11) = -2.*REAL(DPHIDX) 

CP 1 ( 1 i = - 2 . I M AG I DP H I DX ) 

IF { IPRINT ,NE. 1 ) RETURN 
WR ITE ( 6 » 90 ) 

WRITE! 6 T 95) 

WRITE (6, 100) X( U ,0PHIS1 ( 1 ) r DPHIS2( I ) ,PHI( 1 ) , 

1 R( 1 ) , DPHIORt OPHIDX 
WRITE<6jllO) CPR( LlfCPKl ) 

RETURN 

90 FORMAT (' 1 T 1 0 X ‘ , T2 2 D ( PH I )/ D ( S ))% T 46 , 

1 *0( PHU/DIS2) * ,T70t'PHI 5 ) 

95 FORMAT { 5 « ,T10, * R» ,T23, *D(PHI ) /DR* ,T47, * D{ PHI l/DX* ,/ ) 

100 FORMAT ( 3H I , 2X , F 8 .4 7 4X , 3 ( F 1 2 . 5 , F9 . 5 , 3X ) , / * ' ,4X,F3.4, 

I 4X,3(F12.5,F9.5,3X) ) 

110 FORMAT ( ' 0 1 » 10X » 8 INITIAL VALUE: INNER SURFACE OF', 

1 • OUTER CYLINDER* , /* * »30X » * C P( RE AL) = *»F8.4,10X, 

2 'CP(IMAG) = *,F8.4t//) 

END 

THIS IS SUBROUTINE COMPXR 
SUBROUTINE COMPXR 

SUBROUTINE COMPXR COMPUTES THE <X,R) COORDINATES Or A 
GRID POINT. 

DIMENSION X ( 400 ) , R ( 400 ) , XPT ( 400 ) , CPR (400 ), CP I (400) 
COMMON/ BLCK1 / FSTRMN* DELTAS ♦ DX» DH » DSTST R , RO,RI , RE DFRQ 
COMMON/ BLC KB / NGRDFN , NPTS » LCOUNT » IHAVEP ,N,M ,M1 , 1 PR INT 
COMMON/ BLCK5/ X , R , XPT , CPR » CP I 

I = IHAVEP + 1 
J = IHAVEP 

IF THIS IS THE FIRST POINT ON A NEW RIGHT RUNNING MACH 
LINE, GO TO 1. 

IF ( J .EQ. 0 ) GO TO 1 

X ( I ) = X ( J ) + DX 
R(I) = R(J) -DH 
RETURN 

1 R ( I ) = RO 

X( I) = X( 1) + DSTSTR 
RETURN 
END 

C THIS IS SUBROUTINE MACHLN 
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SUBROUTINE MACHLN 



SUBROUTINE MACHLN COMPUTES THE VALUES 0^ THE FLOW 
, H r£k D .,R!^ T l T l ES ALONG THE INITIAL RIGHT-RUNNING MACH 
LINE WHERE PHI AND DPHIS2 ARE ZERO. 



COMPLEX AyC 

COMPLEX DPHIDR, DPHIDX 

COMPLEX PHI (400) , DPHI SI ( 400 ( ,DPHIS2( 400 I 

DIME NS I ON X ( 400 I ,R(400! ,XPT(400 I < CPR <400 ) » CP I (400 ) 

COMMON/ BLCKI / F S T R M N , D E L T A S , D X y D H , D S T S T R , R 0 y R I , R E D F R C 

COMMON/ BLCK2/ A I , A 2 , A3 , A4, RN» R M , RM1 , F ACT I , F ACT2 

COMMON/ 6LCK3/ NGRDFN, NPTS r L COUNT , IHAVEP ,N, M, Ml, I PRINT 

COMMON/ BLCK4 / PHI , DPHI SI ,DPHIS2 

COMMON/ BLCK5/ X , R , XPT , CPR , CP I 



I = IHAVEP 1 
J = IHAVEP 



CALCULATE THE COEFFICIENTS ASSOCIATED WITH THE RIGHT- 
RUNNING MACH LINE. 

LINE. 



A A 1 = 1 . - A 1 / R ( I ) 

A = C. MPLX ( AA1 , A3 ) 

Cl = 1. ^ AI/RU ) 

C = DPHISI ( J)*CMPLX(C1 s-A3J 

PHI (I) = (O.yO.) 

DPHI S2 ( I ) = {0. ,0. 5 
DPHISK I I = C/A 

DPHIDR = A4*0PHI SI ( I ) 

DPHI DX = A2*DPHI Si ( I ) 

IF ( IPRINT .NE. 1 ) RETURN 

WRITE (6, 100) X( I) , DPHISK II,DPHIS2(n,PHI(I), 

1 R( I ) , DPHIDR, DPHIDX 
WRITE! 6, 1 1 1 ) 

RETURN 

100 FORMAT ( 3H M , 2X , F 8 .4 , 4X , 3 ( F 1 2 . 5 , F9 . 5 , 3X I, / « *,4X,F8.4, 
1 4X»3(F12.5 f F9.5 t 3X) ) 

1 1 1 FORMAT ( '0* ) 

END 

THIS IS SUBROUTINE RADI 

THIS IS SUBROUTINE RADI 
SUBROUTINE RADI 

SUBROUTINE RADI COMPUTES THE FLOW FIELD QUANTITIES AT 
THE SURFACE OF THE INNER RADIUS WHERE THE BOUNDARY 
CONDITION IS DEPENDENT ON N , THE CIRCUMFERENTIAL 
MODE NUMBER. IF N IS: 

EVEN — DPHIDR = 0, AND THUS DPHISI = DPHIS2. 

ODD — CP = 0. 

COMPLEX A ,B,C, TANR,DPHIDR, DPHI DX , DELT A , G AMM A * I K 
COMPLEX PHI (400) , DPHISK 400) ,DPH IS 2(400) 

D I ME NS 1 ON X ( 400 ) , R ( 400 ) » XPT ( 40 0 ) , CPR (s-00 ) , CPI (400) 
COMMON/ BLCKI/ F S T RMN, DELI AS , DX , Dri , DS TSTR , RQ , RI , REOFRC 
COMMON/ BLCK2/ A1 ♦ A 2 , A3 , A4 , RN , R M , RM 1 , F ACT 1 1 F ACT2 
COMMON/ 8LCK2 / NGRDFN, NPTS, LCOUNT , IHAVEP ,N,M,M1, IPRINT 
COMMON/ 8LCK4 / PHI , DPHISI ,DPHIS2 
C0MM0N/BLCK5/ X, R , XPT , CPR , CP I 

I = IHAVEP * I 
J = IHAVEP 

CALCULATE THE COEFFICIENTS ASSOCIATED WITH THE RIGHT- 
RUNNING MACH LINE. 
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THE FOLLOWING VARIABLE, IK, WAS FOUND TO BE 
UNDEFINED IN THE ORIGINAL REFERENCE 12 PROGRAM, AND • 
WAS THE REASON FOR A NUMERICAL INSTABILITY AFTER THE 
FIRST REFLECTION. 

IK=CMPLX< 0 .0, REDFRQ) 

AA1 = 1 . - Al/R{ I ) 

A = CMP LX ( AA1 A3) 

B 1 = A1/RII5 - 0.25*DELTAS*DELTAS*(FACT1 - 
I FACT2/(R(I)*R(I»n 

e = CMP LX <B 1, A3 ) 

Cl = 1. <- A 1 / R ( J ) 

C 2 = — A 1 / R { J ) 4- 0. 25*DELTAS*DELTAS*< FACT 1 - 
1 FACT2/I R< I)*R< I) ) ) 

C3 = ( 2 . * F AC T 1 - FACT2/(R< J } *R ( J)| - 
1 FACT 2/ { R ( I)*R(I) ) > *DELTAS*0 . 5 

C = DPHIS1 <J>*CMPLX(C1,-A3) 4- DP HI S2 ( J » *CMPLX( C2 , - A3 ) 

1 + PH 1 ( J ) * CM PL X ( C3 » 0 . ) 

IF ( MOD (N,?} .EQ. 0 ) GO TO 1 

APPLY THE BOUNDARY CONDITION FOR ODD VALUES OF N. 

CP = (IK) PHI 4- DPHIDX = 0. USE THIS B.C., THE FINITE 
DIFFERENCE EQN. FOR DPHIDX, AND THE EQN. FOR THE RIGHT 
RUNNING MACHLINE TO SOLVE FOR DPHIS1 AND DPHIS2. 

DELTA = A 2 + 0.5*1 K*DELTAS 

GAMMA = IK-M PHI C J ) 4- 0 . 5*0 ELT AS*DPH I S2 ( J ) } 

SOLVING FOR DPMI Si AND DPHIS2 

DPHIS1II) = ( C 4- B*GA MMA/DELT A )/( A - B*A2/DELT A ) 
DPF I S2 ( I J = ( C - A*DPHIS1(I) ? / B 

DPHIDX = A 2* ( D PH I SKI) 4- DPHIS2U) } 

D PH I DR = A 4* ( DPHISKU - DPHlS2(Ii ) 

GC TO 2 

APPLY THE BOUNDARY CONDITION FOR EVEN VALUES OF N. 

1 DPHIS1 ( I > = C/( A4-B) 

DPHIS2 ( 1 } = DPHIS1 ( I ) 

DPHIDR = (0.,0.) 

CPHIDX = A 2* ( DPH I S 1 ( I ) 4- DP H I S 2 ( I ) I 
PHI . 

2 PHI ( I ) = PH I ( J ) 4- 0 . 5* I DPH I S2 ( J ) + DPHIS2II ) ) *DELT AS 

IF ( IPRINT .NE. 1 1 RETURN 

WRITE (6, 100) XU) ,DPH IS1( I) »DPHIS2( I ) ♦ PHI l I ) , 

1 R ( I ) , DPHIDR, DPHIDX 
RETURN 

100 FORMAT ( 3H RI , 2X , F8 . 4 , 4X , 3 ( F 1 2 . 5 , F9 . 5 , 3X ) , / ‘ ‘,4X,F8.4, 

1 4X,3 ( F12.5,F9. 5,3X) ) 

END 

THIS IS SUBROUTINE RAD2 



SUBROUTINE RAD2 

SUBROUTINE RAD2 COMPUTES THE FLOW FIELD QUANTITIES AT 
THE INNER SURFACE OF THE OUTER CYLINDER, WHERE THE 
FLOW TANGEMCY CONDITION PRESCRIBES DPHIDR. 



CCMPLEX D , E, F, TAN R 
COMPLEX DPH I DR, DPHIDX 

COMPLEX PHI (400) ,DPHIS1(400) ,DPHIS2(400) 

DIMENSION X ( 400 ) , R ( 400 ) , XP T ( 400 ) ,CPR(400) , CPI (400) 
COMMON/ BLCKl / FSTRMN, DELTAS, DX , DH , DSTSTR ,R0 ,R I , RE DFRQ 
CCMMCN/ 3LCK2/ A1 , A2 , A3 , A4 , RN, RM , RM1 , FACT I , F ACT2 
COMMON/ BLCK3/ NGRDFN , NPTS , LC.OUNT , I HA VEP , N , M , MI , I P R I NT 
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COMMON/ BLCK4/ PH I , DPH I SI T DPHI S2 
COMMON/ BLCK5/ X , R » XP T , CPR , C P I 
OATA PI/3. 141593/ 



I = IHAVEP 1 
J = IHAVEP + 2 
K = LCOUNT 
ARG=RM*PI*X< I) 



CALCULATE THE COEFFICIENTS ASSOCIATED WITH THE LEFT 
RUNNING MACH LINE. 

0 1 = A1/R(I) <- 0.25*DELTAS*DELTAS*(FACTl - FACT2/(R(I) 
1 * R ( I } ) ) 

D = CMP LX { D1 * -A3 J 
El = -1 . - A 1 / R ( I ) 

E = CMPLXIEl t — A 3 J 

FI = -Al/RU) - 0.25*D£LTAS*DELTAS*( FACT1 - FACT2/ 

1 CR!I)#R(im 
F 2 = -1 . 4 - A 1 / R ( J ) 

F3 = <2.*FACT1 - F ACT 2/ ( R ( I ) *R ( I ) ) - FACT2/(R!J)* 
l R( J) })*DELTA$*0.5 

F = DPHISH J)*CMPLX(F1,A3) + DPHI S2( J1 *C MPLXt F2, A3) - 
1 PHI ( J)*CMPLX( F3,0. J 

CALCULATE THE FLOW TANGENCY CONDITION - T ANR 
BOUNDARY CONDITIONS FOR PANEL FLUTTER 



THE FOLLOWING THREE CARDS WERE MODIFIED FROM THE 
ORIGINAL REFERENCE 12 PROGRAM . 



ALPHA-. L 
T ANR1 = ALPHA 

TANR2=REDFRQ*X(I )*ALPHA 
T ANR = CMPLX ( I ANRI * TANR2 J 

DFHJSl(I) = ( F TANR*E/A4 ) / ( 0 <- E ) 

DPFlS2Ui = ( F - TANR*D/A4 )/< D 4- E ) 

DP H I DR = A4* ( D°H I S 1 ( I ) - DPHIS2!I)> 

DPHI DX = A2*(DPHIS1( 1) «• DPHIS2U)) 

INTEGRATE ALONG THE LEFT RUNNING MACH LINE TO FIND PHI 

PHKI) = PHI { J ) 4- 0 .5*1 DPMI SI ( J I + DPHISK I ) ) ^DELTAS 

CALCULATE THE PRESSURE COEFFICIENT AT THIS GRID PG I NT. 



X FT ( K ) 
CPRIKJ 
CPHKI 



-2.#<-REDFRQ*AIMAG(PHI< I) ) + REAL! CPH IDX ) ) 
-2.*( REDFRQ*REAL < PH III) } 4- AI MA G! D PH IDX )) 



IF ( IPRINT .ME. 1 ) RETURN 

WRITE! 6, 1.00) XII) , DPHI SI ( I) , D PH I S2 ( I J » PH I! I ) , 

1 R(I) ,DPHIDR.DPHIDX*TANR 
WRITE(6,I10} CPR(K) » C PI IK) 

RETURN 

100 FORMAT ( 3 FI RO, 2X t F8 .4r 4X t 3< FI 2 . 5 » F9.5 « 3X ) « / 1 *,4X,F8.4, 

I 4X,3(F12.5»F9.5f3X) ) . 

110 FORMAT! < 0 « ,10X,' INNER SURFACE OF OUTER CYLINDER' ,/' 

1 30X, 'CP(REAL) = • ,F8.4,10X, 5 CP( I MAGI = *»F8.4,//) 

END 

THIS IS SUBROUTINE GENFPT 



SLB ROUT I NE GENFPT! IER) 

SUBROUTINE GENFPT COMPUTES THE FLOW FIELD QUANTITIES 
AT A GENERAL FLOW FIELD POINT. 



COMPLEX PHIL 
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COMPLEX A,B»C,D,E, F, DENOM 
COMPLEX D PHI OP » DPH I DX 

CCMPLEX PHI! 400 ) ,DPHISU400) ,DPHIS2(400) 

DIMENSION X( 4001 ,R (400) ♦ XPT { 400 ) , CPR < 40 0 ) , C PI ( 400 ) 
COMMON/ BLCKl / FSTRMN, DELTAS ,DX, DH,0STSTR, R0, RI , RE DERG 
COMMON/ BLCK2/ A i , A 2 , A 3 , A4,RN, RM, RM1 , F ACT 1 , F ACT2 
COMMON/ BLC K3 / NGRDFN , NPTS , L COUNT , IHAVEP , N, M , ML , I PRINT 
COMMON/ BLCK4/ PHI ,DPHISl ,DPHIS2 
COMMON/ B LCK5 / X , R , XPT , CPR , CP I 

I = IHAVEP 4- 1 
J = IHAVEP 
K = IHAVEP *■ 2 

COMPUTE THE COEFFICIENTS QF DPHISI AND DPHIS2. 

FACTS = FACT?/! R( I)*R{ I ) I 
FACT 4 = FACT2/(R(K)*R(KI) 

FACT5 = FACT2/I R( J >*R( J) ) 

A , B AND C ARE ASSOCIATED WITH THE RIGHT RUNNING MACH 
LINE, WHILE D, E AND F ARE ASSOCIATED WITH THE LEFT. 

AA1 = 1 . - A I / R ( I ) 

A = CMPLX ( AA1 , A3 ) 

81 = -A A I 4- 1. - 0 .25*DcLTAS*OELTAS*< FACT1 - FACTS) 

B = CMP LX { 81, A3 > 

Ci = 1. <- A 1 / R { J ) 

C2 = -Cl 4-1.4- 0 . 2 5*DELTAS*DFLTAS*( FACT 1 - FACT 3 ) 

C3 = ( 2 . * FACT 1 - FACT 3 - FACT5 ) *DELTAS*0 .3 
C = DPHISI! J)*CMPLX!Ci,-A3) * DPMI S2{ J) *CMPLX(C2 , -A3 ) 
1 4- PHI ( J1*CMPLX(C3,0. j 

D1 = Al/RU) 4 - 0 .25*DELTAS*DELTAS*( FACT 1 - FACT3 ) 

D = CMPLXtDl ,-A3 ) 

El = ~i . ~ Al/R( I ) 

E = CMPLX! El , - A3 ) 

FI = —A 1 / R ( K ) - 0. 2 5* DELTAS* DELTAS*! FACT 1 - FACT 3 ) 

F 2 = -I. 4- A 1/R ( K ) 

F2 = ( 2 . * F AC T 1 - FACT3 - F ACT4 ) *DE IT AS*0 . 5 
F = DPHISl(K)*CMPLX(Fl» A3) + D PH I S2 ( K ) *C MPLX( F2 ? A 3 ) - 
1 PHI (KJ *CMPLX( F3,0. J 

USE CRAMERS RULE TO SOLVE FOR DPHISI AND DPHIS2 
PROVIDED THAT THE DETERMINANT OF COEFFICIENTS IS NON- 
SINGULAR. 



DENOM = A*E - D*B 

IF ! CABS! DENOM ) ,LT« l.E-05 ) GO TO I 
DPHISI! I) = ( C*E - B*F) /DENOM 
DPHIS21I) = ( A*F - C*D ) /DENOM 
DFHIOR = A4*< DPHISI (I) - DPHIS2(I)) 

DPHIDX = A2* (DPHISI! I ) 4 - DP H IS 2 ! I ) ) 

INTEGRATE ALONG THE RIGHT RUNNING MACH LINE TO FIND 
PHI . 



C 



C 



PHI ! I ) = PHI(J) 4- 0. 5*( DPHI S2 ( I ) 4 - DPHI S 2 ( J ) ) *DE LTAS 
PHIL = PHI(K) 4- 0.5* (DPHISI (K) + DPH I S 1 ( I )) *UELT A S 

IF ( IPRINT .NiE. 1 ) RETURN 

WRITE! 6, 100) X! I ), DPHISI! I ) ,DPHIS2!I),PHI( I), 

1 R( I ) ,DPHIDR, DPHIDX, PHIL 
WRITE! 6, ill) 

RETURN 



1 IER = 1 

WRITE(6»110) DENOM 
RETURN 

100 FORMAT ( 3H G , 2X , F 8 .4, 4X , 3 ( F 1 2. 5 , F9 . 5 ,3X ) , / ' 

I 4X,3 (FI 2. 5,F9.5,3X» ) 

110 FORMAT ( * 0 * ,///,* * , 35X , ' ABNORMAL TERMINATION: 



1 



•DETERMINANT OF COEFFICIENTS IS SINGULAR', 



, 4X , F 8 . 4 , 

• , 

/' 0« , 50 X, 
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2 'DENOMINATOR 
111 FORMAT! ' O') 

END 
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THIS IS SUBROUTINE QSF 
SUBROUTINE QSF 
SUBROUTINE QSF ( H » Y , Z , ND IM ) 

CALL QSF <H,Y,Z,NDIM) 

DESCRIPTION OF PARAMETERS 

H - THE INCREMENT OF ARGUMENT VALUES. 

Y - THE INPUT VECTOR OF FUNCTION VALUES. 

Z - THE RESULTING VECTOR OF INTEGRAL 

VALUES, l MAY BE IDENTICAL WITH Y. 
IDENTICAL WITH Y. 

NDIM - THE DIMENSION OF VECTORS Y AND Z. 



REMARKS 

NO ACTION IN CASE NDIM LESS THAN 3. 



DIMENSION Y( 1 ) t Z ( 1 ) 

HT = .333333 3* H 
I F ( ND I M- 5 ) 7, 8; 1 

NDIM IS GREATER THAN 5. PREPARATIONS OF INTEGRATION 
LCCP 



1 SUf*l=Y( 2 ) *Y( 2) 

SUM 1 = $U Ml * SUM1 
SUM1=HT*( Y( I ) +SUM1+Y ( 3) ) 

AUX1=Y( 4UY14! 

AUX1=AUX1 -s-AlJXI 

A UX 1 = S U M i + HT * ( Y < 3 ) *AUX1+Y< 5 ) ) 

AUX2=HT*( Y( I ) t-3. 875* ( Y( 2 )+Y( 51) +2,625*( Y (3 > + Y(4) ) 4 - 
1 Y ( 6 i ) 

SUM2=Y < 5 ) +Y( 5 ) 

SU^2= SUM2 + SUM2 

SUM2- AUX2-HT * ( Y ( A ) * SUM2 4-Y ( 6 ) ) 

zm=o. 

AUX= Y l 3 ) + Y ( 3 ) 

ALX=AUX+AUX 

Z(2»=SUM2-HT*( Y( 2H-AUX+Y(4) ) 

Z ( 3 ) = S U M I 
Z ( 4 ) = S U M 2 
I F ( ND I M-6 ) 5 , 5 , 2 

INTEGRATION LOOP 

2 DO 4 1 = 7, NDIM, 2 
SUM1=AUX I 
SUM2=AUX2 

AUX l-Y ( 1-1MY! I-l ) 

AUX1 = AL'Xl + AUX1 

AUX1=SUM1<-HT*(Y( 1-2 ) *- AUXi+Y ( I ) ) 

Z( 1-2 ) = SUM 1 
I F < I- N D I M ) 3 , 6 , 6 

3 AUX2 = Y { I ) +■ Y ( I ) 

AUX2=A UX2 <-AUX2 

AUX2=SUM2«-HT*( Y( I-l ) +AUX2*Y< I + l ) ) 

4 U I-l ) = SUM2 

5 Z ( NDI M-l ) = AUX1 
Z ( ND I M ) = AUX2 
RETURN 

6 Z( NDI M-l )=SUM2 
Z < N D I M ) = A U X I 
RETURN 

END OF INTEGRATION LOOP 

7 I F ( NDIM-3 ) 12, lit 8 

NDIM IS EQUAL TO 4 OR 5 

8 SUM2=1.125*HT*(Y<1 ) +Y( 2 ) +Y l 2 ) *Y< 2 H-Y ( 3 ) + Y( 3 KY( 3 ) + 



8 4 



oo 



1 Y ( 4 ) ) 

SUP 1=Y ( 2 I +Y ( 2 ) 

SUP 1 = SUM1 + SUHl 
SUM1=HT* { Y ( 1 )*SUM1+Y< 3) ) 

ni )=o. 

AUX 1 = Y { 3) «• Y( 3 ) 

AUX1=AUXH-AUX1 

Z< 2)=SUP2-HT*( Y( 2 ) *AUXH-Y( 4) ) 

I F ( ND I M-5 ) 10,9,9 
9 A U X 1 = Y ( 4 j « Y ( 4 ) 

A U X I = A UX 1 -?■ A UX 1 

Z { 5 ) = S U M i +• H T * < Y ( 3 ) *- A U X l +• Y ( 5 } ) 

10 Z(3)=SUMi 
Z ( 4) = SUM2 
RETURN 

ND I M IS EQUAL TO 3 

11 SUM1=HT*{ ] . 2 5*Y ( 1 H-Y( 2 ) +Y ( 2 i- . 25*Y (3) \ 
SUM2=Y { 2 > t- Y ( 2) 

S U M 2= S U M 2 + S U M2 
l { 3 ) =HT * (Y(l) +SUM2+Y ( 3 ) } 

Z ( 1)=0. 

Z( 2)=SUMl 

12 RETURN 
END 
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The outer flowfield program, from reference 13, consists 
of 5 subroutines and a main program to manipulate the 
subroutines. The action of each particular subroutine is 
indicated in the program listing. Inputs required to the 
program must be on two cards, and are as follov.'s: 

DATE 



MACH Mach number. 

RADIUS Radius to .length ratio. 

K Reduced frequency. 

MFREQ Axial mode number. 

NR The r in Qmr, the generalized aerodynamic force. 

NFREQ Circumferential mode number. 

FINGSD Grid fineness, essentially the number of grid 
points along the cylinder surface. 

The run date is typed onto the first card. The second card 
contains the above parameters in the following format: 
F20.8, 2R10.3, 415. Df the above parameters, in the 

modified program, MPRFQ, NR, end NFREQ must be set to 1, as 
the case under consideration is not panel flutter but a 
cylinder at angle of attack. Other parameters are 

arbitrary. 
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